Detection of partial polarization in Compton 

scattered photons 
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Abstract 



It has been recently proposed [1] to use polarization of Compton scattered 7- 
rays to improve the imaging performance of Compton telescopes. Building upon 
that work, we detected the aforementioned polarization in a sample of 1.836 MeV 
7-rays from the LXeGRIT Compton telescope. Here we present the measurement, 
together with detector oriented considerations on the application of the principle to 
a realistic Compton telescope. 
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Introduction 



In a recent work [1] it has been shown how polarization of Compton scattered 
7-rays retains information about the source location, and it has been sug- 
gested that it provides a very significant improvement in the reconstruction 
(imaging) of a 7-ray source, the most common imaging technique being based 
on Maximum Likelihood algorithms (see for example Ref. [2]). 
The basic idea can be outlined as follows: given an unpolarized 7-ray source, 
let's assume that two Compton scatters followed by full absorption are de- 
tected. The first scattered photon will be partially polarized in a direction 
perpendicular to the scatter plane; the second scattered photon will be prefer- 
entially in a plane perpendicular to the polarization direction, i.e. in the plane 
containing the source, the first interaction site and the second interaction site. 
It is clear that a fine grained detector is needed in order to be able to detect 
each interaction. LXeGRIT 1 is a good example of such a detector and polar- 
ization has been detected in a sample of 3-site events from 1.836 MeV 7-rays. 
The measurement is described in detail in Sec. 1. We did not try to assess 
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directly the impact of polarization in constraining the source position, since 
it would require a major effort to be implemented in the current LXeGRIT 
imaging procedure. We try, anyway, to clarify what is needed to make this 
approach effective in Sec. 2. 

When coming down to realistic detectors, the field of Compton telescopes is 
nowadays rather fluid, without much agreement within the community about 
the detector of choice. The sheer lack of such a canon makes it difficult to bring 
even bright ideas from the realm of wishful thinking to practical realization. 
To give an example, it was proposed [4] to track the recoil electron in the first 
Compton scatter, which is extremely troublesome since tracking ~MeV elec- 
trons requires a low density, low Z active medium, clearly conflicting with the 
requirement of stopping power needed for 7-rays. A realistic design to realize 
such an idea is not available to date. 



1 Measurement 

For a concise and clear exposition of the general problem, we refer to [1]. 
LXeGRIT is a Compton telescope which exploits a liquid xenon time projec- 
tion chamber (LXe TPC) to detect MeV 7-rays, imaging each interaction in 
the fiducial volume with an accuracy better than 1 mm on each coordinate; 
quite naturally the TPC provides a built-in Cartesian reference frame. Data 
from an 88 Y 7-source, which emits photons at 0.898 and 1.836 MeV, have 
been used for this measurement. The source was placed on top of the TPC, 
at a distance of 2 m along the z axis and centered at x=0 and y—0 in the 
TPC reference frame. The 1.836 MeV line was selected for it provided a better 
detection efficiency. Thus we experimentally deal with 

1. direction of 7-rays coming from the source, which, in the TPC reference 
frame is (0,0,1); 

2. Ei, Xi, i/i, Zi for the i th interaction, z=l,2,3; 

3. E 1 +E 2 +E 3 =1.836 MeV, because we select events in the full energy peak. 

For each event, the three interactions have been sequenced using Compton 
kinematics [5] and only events which give the correct source location are then 
selected; in this way the contamination due to incorrectly sequenced events is 
reduced to a truly negligible level. 

In order to study the polarization of the scattered 7-rays, we change to the 
reference frame described in [1], i.e. the z axis in the direction from the second 
interaction to the first one and the y — z plane defined by the three interac- 
tion sites. In the new reference frame the vector from the first to the second 
interaction is V12 = (0,0,-1) and the vector from the second to the third 
interaction is V23 = (0, sin #123, cos #123); #123 is the second scatter angle. The 
x axis of the new reference frame is now given in the TPC reference frame as 
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the cross product 



x = v u ' x v 23 ' 



properly normalized; primed vectors are in the TPC reference frame. 

The source position, in the new reference frame, is defined as (Eq. 5 in [1]) 

r = (sin #012 cos 0, sin 6> i 2 sin 0, cos #012) (2) 



where 6*012 is the first scatter angle and is the azimuthal angle in the x — y 
plane; is the sensitive variable in this measurement. The x component of 
r is found projecting the source direction (0,0,1) on the x axis as defined in 
Eq. 1, and the y component then follows. Once r 0x has been determined, cos 
is given and we can look for an excess in the distribution, favoring the plane 
containing the source, the first interaction site and the second interaction 
site (012-plane). Here we shift of 90°, so that the 012-plane corresponds to 
0=90°, while in [1] it corresponds to 0=0°, 180°. As shown in [1], this excess in 
the distribution is expected to be small for 1.836 MeV 7-rays; as expected 
in polarization phenomena (and confirmed in Fig. 2 in [1]), it is supposed to 
be maximum for # 12 3 ~ 90° and negligible for # 12 3 < 30° and # 12 3 > 150°. 
To emphasize the excess in the distribution, we divide our sample in two 
sub-samples: one with 66° < #123 < 114° and the other one with 6>i 23 > 
114° or 9i23 < 66° (the precise choice of these numbers is not critical). Our 
expectation is now to detect polarization in the first sub-sample while the level 
of polarization should be negligible in the second one. The distribution for 
the two sub-samples is shown in Fig. 1-left. We define a modulation as 

modulation = — -2- (3) 

01 + 02 



where is the distribution binned, in this case, in intervals of 10°; the 
subscripts indicate the sub-sample. The modulation is shown in Fig. 1-right, 
error bars from statistical errors; a significant excess is visible in the region 
= 90° ± 30°, as expected. 

We can look at the problem in a slightly different way, taking Eq. 10 in [1] as 
our starting point; it can be rewritten as 

(3 + P~ 1 -2 sin 2 ol2 (/( - (I[ - IQ cos 2 0) 
9oi2) = MPT^ -2sin 2 6W) (4) 
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Fig. 1. Left: 4> distributions for the two sub-samples 66° < #123 < 114° (continuous 
line) and #123 > 114° or #123 < 66° (dashed line). Right: modulation vs. 4>. 

where (3 is the ratio of the photon energy after the scatter to its energy before 
scattering, i.e. 

1 0.511 + E 7 (l - cosO Q12 ) 
E 1 = 1.836 MeV in the present case; 

p & + P'- 1 



2(/3' + (3'~ 1 -sm 2 e 123 ) 

_ P' + P'- l -2sm 2 9 123 
2 2(/3> + /3>- 1 -sm 2 6 123 ) 
where primed quantities refer to the scattered 7-ray and E^—E\ + E 2 in /?'; 
here (ft is defined as in [1]. 

For unpolarized photons I[ — I 2 — 0.5 and Eq. 4 is reduced to 

~(a \ £ ± P' 1 - sin2 9 ° 12 (k\ 
P{0O12) ~ MPTF 1 -2sin^ 012 ) (5) 



which does not depend on 0, as expected. We can therefore compare p(<p, #012) 
to p(0oi 2 ) and look for some discrepancy. Again, it is advisable to split the data 
in two sub-samples and pick up the differences between them. The results for 
the aforementioned sub-samples are shown in Figs. 2, 3. Fig. 2-left shows 
K0i ^012) vs - c °s ^012 for the supposedly unpolarized sub-sample (2), Fig. 2- 
right for the supposedly polarized sub-sample (1); ^(^012), obviously the same 
for the two different samples, has been superimposed (continuous line). In an 
attempt to quantify how much discrepancy from the unpolarized case is in the 
two sub-samples, the histograms of the residuals are shown in Fig. 3; the dashed 
line is for sub-sample 2, the continuous line for sub-sample 1. The supposedly 
polarized sub-sample clearly shows a more pronounced discrepancy from the 
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unpolarized case and its residual distribution is shifted toward negative values, 
as expected. 
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Fig. 2. p(cj), 6012) vs. cos 0oi2 compared to ^(#012) vs. 0036*012 (continuous line). Left: 
sub-sample 2, right: sub-sample 1. 
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Fig. 3. p((p, #012) — p(#oi2)j i- e - residuals from Fig. 2; the dashed line is for 
sub-sample 2, the continuous line for sub-sample 1. 

2 Discussion 

We have presented the detection of polarization in a sample of Compton scat- 
tered 7-rays, following the suggestion in [1]. It was done for a relatively high 
energy (1.836 MeV), where the effect is expected to be small, because it offered 
the larger sample. 

To quantify the modulation as defined in [1] is well beyond the scope of this 
work, as well as to quantify its impact on source imaging. A rough but solid 
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estimate can be given answering the following question: which fraction of imag- 
ing events 2 ends up in the polarized sub-sample, as defined earlier? 
In LXeGRIT, a large fraction of imaging events is detected with two inter- 
actions only, and for this 2-site events no polarization of the scattered 7-ray 
is detectable, as obvious; at 1.836 MeV, 2-site events make up ~80% of the 
imaging events. 2-site events include a large fraction of the events with a large 
(~90°or larger) scatter angle, which are most favorable to detect polarization. 
Of the remaining 3-site events, about 50% are in the polarized sub-sample, 
i.e. ~10% of the total. This figure gives a pre-factor that should be included 
when quoting the modulation along the event circle; the ratio 3-site events to 
imaging events is further reduced at lower energies. 

It is clear at this point that polarization of Compton scattered 7-rays is a 
small effect for a detector like LXeGRIT, and this conclusion can be easily ex- 
trapolated to other detectors using dense, high Z active media (e.g. Ge) and a 
mm position accuracy. On the other side, polarization has been detected using 
a detector by no means optimized to do so, possibly indicating that to use it 
to constrain the event circle in Compton imaging is within reach of a realistic 
Compton telescope. 

The main limitation in LXeGRIT does not come neither from the relatively 
high energy threshold (150 keV) or from its mm spatial accuracy (very hard to 
improve in any realistic, sizeable 7-ray detector), but from LXe itself, which 
has a large cross-section for photoelectric absorption due to Zx e = 54, exceed- 
ing the Compton cross-section below 300 keV. LXe has a density of ~3 g/cc 
and an attenuation length well below 1 cm for photon energies of 200 keV or 
less, so that many interactions are too close to be spatially resolved. These 
properties, which allow to build an extremely compact and efficient Compton 
telescope out of LXe, pose a significant limit on the event multiplicity in LXe, 
when a multiplicity > 3 is required for detecting polarization of Compton 
scattered 7-rays. 

Without abandoning the terrain of noble liquids, a sound alternative would be 
liquid argon (LAr), since Z^ r = 18 and its density is ~1.4 g/cc. A LAr TPC 
is a viable alternative to a LXe TPC: LAr TPC's with the fiducial volume 
necessary for a Compton telescope have been built [6], LAr is purified more 
easily than LXe, LAr is relatively inexpensive and, last but not least, LAr 
allows a better energy resolution than LXe [7] . 



Conclusions 

In this work, stimulated by [1], we have presented the detection of polarization 
for Compton scattered 7-rays at 1.836 MeV; data from the LXeGRIT Compton 

2 For an imaging event, full containment and at least two interactions are required. 



6 



telescope have been used. On the one hand, polarization has been detected; 
on the other hand, it is clear to us that such an effect is small when any high 
density, high Z material (LXe in our case) is used as active medium. Based on 
our experience with LXeGRIT, LAr is proposed as a promising alternative. 
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